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Section  1 
INTRODUCTION 

Fluid  Analysis  P^rogram,  _3-Dlmensions  (FLAPS)  is  a three- 
dimensional  finite  element  code  developed  for  the  purpose  of  predicting 
the  combustion  instability  of  solid  propellant  rocket  motors.  The 
complete  code  consists  of  four  phases: 

a)  A pre-processor  phase  (a  three-dimensional  finite  element 
mesh  generator). 

b)  The  acousto-modal  analysis  (calculation  of  the  acoustic 
frequency,  mode  shape,  and  corresponding  velocities). 

c)  The  potential  flow  analysis  (calculation  of  gas  flow 
velocities). 

d)  The  determination  of  mode  stabllity/lnstabillty  (evaluation 
of  stability  integrals) . 

The  second  phase  of  the  code  provides  for  the  coupling  of  the 
acoustic  response  of  the  rocket  cavity  with  that  of  the  solid  propellant 
grain  for  the  purpose  of  obtaining  a more  accurate  value  (than  would  be 
obtained  by  considering  "acoustically  hard"  cavity  boundaries)  of  the 
acoustic  frequency,  and  for  calculating  the  amount  of  structural 
damping  provided  by  the  propellant  grain.  The  solid  propellant  can  be 
modeled  as  a virtually  incompressible  material  having  a frequency 
dependent  complex  modulus.  The  frequency  dependency  of  the  solid 
propellant  modulus  may  require,  in  some  cases,  the  employment  of  an 
iterative  technique  to  obtain  the  coupled  frequency  and  the  rate  of 
structural  damping. 

Because  modeling  of  the  solid  ^..opellant  greatly  increases  the 
size  of  the  problem,  and  therefore  the  cost  of  running  it,  the  user 
has  the  option  of  a cavity-solid  coupled  analysis  or  of  a cavity 
analysis  alone.  The  principle  of  dihedral  symmetry  is  fundamental  in 
the  development  of  the  code  and  requires  that  the  user  generate  a 
finite  element  mesh  only  for  the  geometry  of  the  smallest  repeating 
cavity-propellant  segment,  the  problem  being  basically  one  in 
cylindrical  coordinates. 

FLAP3  was  developed  primarily  for  the  following  purposes: 

a)  To  provide  generality  and  accuracy  in  the  modeling  of 
complex  combustion  chamber  geometries. 

b)  To  provide  a means  of  predicting  the  damping  of  acoustic 
oscillations  by  the  solid  propellant  grain. 

c)  To  provide  an  integrated  program  designed  solely  for  the 
purpose  of  combustion  instability  prediction  with  ease  of  use. 


TJie  main  features  of  KI,AP3  are  given  In  tlie  following  list  and 
will  be  explained  later; 

a)  It  utilizes  a specially  adapted  three-dimensional  finite 
element  input  data  generation  (FI.ESH3  and  FLSH3P)  package. 

b)  It  utilizes  the  principle  of  dihedral  symmetry  which  enables 
a consideration  of  only  the  smallest  repeating  geometrical  segment. 

c)  It  couples  the  response  of  the  gas  cavity  with  that  of  the 
solid  propellant  grain  to  enable  the  calculation  of  the  frequency  of 
the  coupled  system  and  the  acoustic  wave  damping  provided  by  the 
propellant  grain. 

d)  It  provides  for  modeling  the  propellant  grain  as  a nearly 
incompressible  material  (which  differs  from  the  usual  minimum  potential 
energy  formulation). 

e)  It  utilizes  the  principle  of  condensation,  in  which,  in  this 
case,  the  fluid  pressure  degrees-of-freedom  are  designated  '’master"  and 
the  solid  propellant  displacement  and  mean  pressure  parameter  degrees- 
of-freedom  are  designated  "slave."  This  enables  a major  reduction  in 
the  size  of  the  problem;  the  number  of  equations  is  reduced  from  the 
total  number  of  degrees-of-freedom  of  the  coupled  system  to  the  number 
of  fluid  pressure  degrees-of-freedom. 

f)  It  provides  the  option  of  considering  the  response  of  the  gas 
cavity  alone  (which  raodexs  the  cavity  boundaries  as  "acoustically  hard") . 
This  option  might  be  utilized  in  certain  cases  where  a savings  in 
computer  costs  or  storage  is  a dominant  consideration.  In  this  case 

the  previously  described  condensation  routine  obviously  would  not  be 
nvoked . 

g)  It  calculates  the  three-dimensional  potential  flow  field. 

h)  It  evaluates  the  stability  integrals  for  the  calculation  of 
the  net  driving/damping  coefficient  for  each  acoustic  mode. 


Section  2 

PRE  PROCESSOR  - FINITE  ELEMENT  MESH  GENERATOR 


It  Is  necessary  to  develop  a finite  element  mesh  for  only  one 
repeating  segment  (Figure  1)  of  the  total  cavity-solid  propellant 
rocket  geometry.  This  is  true  because  of  the  employment  of  the 
principle  of  dihedral  symmetry  in  FLAP3.  Although  the  three-dimensional 
element  used  in  FIAP3  for  both  the  cavity  region  and  the  solid  propellant 
is  a tetrahedron  [1],  the  mesh  is  that  of  bricks  connected  at  the 
corners.  Each  brick,  or  quasi-hexahedron,  is  comprised  of  five  basic 
tetrahedra  (Figure  2).  The  breakdown  of  the  quasi-hexahedron  into 
tetrahedra  is  performed  internally;  therefore,  the  mesh  need  be  no 
more  complicated  than  that  shown  in  Figure  3. 


SMALLEST 

REPEATING 

SEGMENT 


Figure  1.  Three-dimensional  cavity-solid  propellant  model. 


The  pre-processor  phase,  or  finite  element  mesh  generator,  FLESH3 
(Fluid  Mesh  Generation,  _3-Dimensions) , is  an  efficient  routine  which 
automatically  creates  the  complete  finite  element  mesh,  consisting  of 
both  cavity  and  solid  propellant  regions,  from  a minimal  amount  of 
input.  Each  repeating  segment  is  sectioned  in  the  longitudinal  or 
z-direction,  with  each  section  comprised  of  a number  of  quadrilateral 
parts  which  are  identified  by  a counterclockwise  listing  of  their  part 
boundary  curves.  Part  boundary  curves  may  be  ellipses  as  well  as 
straight  lines;  their  points  of  intersection  are  designated  by  I,J 
indices  as  shown  in  the  User's  Guide  for  FLESH3  (Figure  4)  which  begins 
on  page  14. 


Section  3 
DISCUSSION 


3. 1 Dihedral  Symmetry 

If  a geometrically  defined  body  Is  comprised  of  Identical 
segments  symmetrically  arranged  with  respect  to  an  axis,  the  degrees- 
of-freedom,  for  a finite  element  analysis,  can  be  transformed  Into 
uncoupled  symmetrical  components,  thereby  greatly  reducing  the  number 
of  equations  which  must  be  solved  simultaneously  [2].  A further 
reduction  occurs  If  each  segment  has  a plane  of  reflective  symmetry. 
Dihedral  symmetry  Is  the  term  applied  to  this  latter  condition.  It 
can  be  seen  from  Figure  1 that  a typical  solid  rocket  geometry  meets 
this  requirement;  therefore,  the  principle  of  dihedral  symmetry  can 
be  employed  in  a three-dimensional  (cylindrical  coordinate)  analysis. 
The  application  of  the  principle  to  this  problem  is  explained  In  detail 
in  Reference  3 and  will  not  be  repeated  here,  but  the  resulting 
analysis  will  be  discussed. 

Employing  the  principle  of  dihedral  s)rmmetry,  FIAP3  computes 
three  distinct  types  of  acoustic  harmonics:  the  zero  harmonic,  the  K 

harmonics,  and  the  M/2  harmonic.  The  zero  harmonic  exists  for  all 
cases  and  is  the  only  harmonic  which  does  exist  for  the  case  of  a 
rocket  having  an  annular  cavity.  The  number  of  po  sible  K harmonics 
is  given  by 


K = 1 J 


(1) 


where  the  harmonic  index  J is  given  by 


J = 


(if  M is  an  odd  number) 


(2a) 


or 


J 


(if  M is  an  even  number) 


(2b) 


where  M is  the  number  of  radial  slots  or  lobes.  For  the  case  of  an 
annular  cavity,  where  radial  slots  do  not  exist,  an  arbitrarily  small 
repeating  segment  (subtended  by  a whole  number  of  degrees)  may  be 
selected  for  analysis.  The  M/2  harmonic  exists  only  when  M is  an  even 
number.  Referring  to  the  geometry  of  Figure  1,  the  zero,  first,  and 
second  harmonics,  for  example,  could  be  calculated.  The  latter  two 
are  both  K harmonics. 
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The  longitudinal  modes  associated  with  each  harmonic  are  calculated 
internally  by  FLAP3,  as  requested  by  the  user.  The  result  of  this 
operation  (acousto-modal  elgensolution)  is  the  natural  circular  frequency 
associated  with  each  acoustic  mode  and  the  corresponding  acoustic 
pressure  distribution  (normalized  acoustic  pressure  at  each  finite 
element  nodal  point)  for  the  smallest  repeating  segment.  The  pressure 
distribution  in  all  of  the  other  segments  is  then  simply  calculated 
automatically  through  the  dihedral  transformation.  The  acoustic 
velocity  components  (constant  for  the  region  occupied  by  each  cavity 
tetrahedral  element)  are  computed  from  the  acoustic  pressure  nodal 
point  values  for  those  cavity  elements  which  are  adjacent  to  the 
solid  propellant  grain  (characterized  in  FLAP3  as  fluid  interface 
tetrahedrons) . 

The  formulation  of  the  complete  three-dimensional  finite  element 
acousto-modal  analysis,  in  which  the  natural  circular  frequency,  the 
acoustic  pressure  distribution,  and  the  element  acoustic  velocities 
are  calculated,  is  given  in  Reference  4 and  will  not  be  repeated  here. 

The  theoretical  finite  element  formulation  which  was  used  in  the 
development  of  FLAP3  is  also  presented  in  Reference  5. 

3. 2 Coupled  Response 

The  presence  of  the  solid  propellant  grain  can  significantly 
shift  the  acoustic  system  frequency  from  that  of  the  gas  phase  alone, 
a portion  of  the  acoustic  energy  being  dissipated  by  the  deformable 
solid  material.  This  effect  can  be  one  of  the  more  significant  sinks 
for  acoustic  energy  in  both  large  and  small  rocket  motors,  the  amount 
of  damping  depending  on  the  grain  geometry  and  mechanical  properties 
and  on  the  acoustic  mode  shape  and  natural  frequency. 

To  evaluate  the  coupled  cavity-solid  propellant  grain  response, 
it  is  necessary  to  model  both  cavity  and  grain  by  the  finite  element 
method.  This  greatly  increases  the  size  of  the  problem  to  be  solved 
from  the  standpoint  of  number  of  initial  degrees-of-freedom.  The 
coupled  finite  element  formulation  [5]  which  is  programmed  in  FLAF3 
is  expressed  in  matrix  form  as: 


T I U ' 

I 

o’ 

0 1-1+J:'. 
10(0 


0 


(3) 


where  [F]  is  the  fluid  inertia  matrix,  [T]  is  the  fluid  compressibility 
matrix,  [K]  is  the  solid  stiffness  matrix,  [M]  is  the  solid  consistent 
mass  matrix,  [U]  is  the  matrix  which  couples  acoustic  pressure  degrees- 
of-freedom  to  solid  displacement  degrees-of-freedom,  {p}  is  the  acoustic 
' pressure  vector,  {A}  is  the  solid  displacement  vector,  (H)  is  the  mean 

2 

pressure  parameter  vector  (to  be  explained  later),  and  X is  the  eigen- 
value of  the  coupled  system. 
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The  srrmrural  damping  can  be  attributed  to  t lie  out-of-pha.se 
response  of  the  solid  propellant  grain  which  is  mi>asured  in  terms  of 
the  complex  shearing  modulus  of  the  grain,  which,  in  turn,  results  in 
a complex  eigenvalue  for  the  coupled  system.  The  imaginary  part  of 
the  complex  eigenvalue  obtained  from  the  e igensoj ut i on  is  the  natural 
circular  frec|uency  of  the  coupled  system  while  the  real  part  is  the 
structural  damping  rate. 

Because  the  ciimplex  shearing  modulus  is  frequency  dependent,  a 
series  of  iterations  may  be  necessary  before  the  accurate  value  of 
complex  modulus  for  input  into  the  program  is  determined. 

3.  3 Propellant  Grain  Modeling 

Because  the  propellant  grain  is  accurately  modeled  only  as 
a nearly  incompressible  material,  the  well-known  standard  Navier  dis- 
placement formulation  would  lead  to  Inaccuracies  in  the  finite  element 
modeling  of  the  grain.  To  avoid  this  situation,  the  solid  finite 
element  formulation  utilized  in  FLAP3  is  that  of  a linear  displacement 
linear  mean  pressure  tetrahedron  [6].  It  is  similar  to  the  Herrmann 
variational  formulation  [7]  which  employs  a linear  displacement 
function  and  a constant  mean  pressure  function.  The  finite  element 
modeling  of  the  propellant  grain  used  in  FLAP3  is  outlined  in  detail 
in  Reference  5 and  will  not  be  repeated  here. 

3. 4 Eigenvalue  Economizer  — Condensation 

The  extraction  of  eigenvalues  and  eigenvectors  is  a far 
more  expensive  operation  than  is  the  solution  of  simultaneous  linear 
equations.  It  requires  roughly  twice  as  long  to  extract  a single 
eigenvalue  as  it  does  to  do  a single  static  analysis.  To  reduce  or 
condense  the  number  of  degrees-of-f reedom  in  the  eigensolution,  the 
following  technique  is  utilized  in  FLAP3.  Further  details  of  the 
method  are  found  in  Reference  8. 

The  original  formulation  of  the  coupled  system  is  given  by 
Equation  (3),  where  the  number  of  degrees-of-f reedom  is  equal  to  the 
number  of  cavity  nodal  point  pressures  plus  the  number  of  solid 
propellant  nodal  point  displacement  components  plus  the  number  of 
solid  propellant  nodal  point  mean  pressure  parameter  values  (one  at 
each  propellant  node)  for  the  analyzed  repeating  segment.  The 
condensed  formulation  is  given  by: 
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where  I is  the  identity  matrix.  The  relationship  between  initial  and 
reduced  degrees-of-freedom  is  given  by: 


Although  one  of  the  most  important  features  of  FLAP3  is  the 
coupled  cavity-solid  propellant  response  analysis  capability,  the 
option  of  a cavity  analysis  alone  is  available  to  the  user.  In  this 
case,  the  cavity-solid  propellant  Interface  is  modeled  as  an 
"acoustically  hard"  boundary.  The  cavity  only  option  would  greatly 
reduce  the  computer  storage  requirements  for  FLAP3  and  would  in  certain 
cases,  suffice.  The  user  need  make  only  two  simple  modifications  in 
the  input  data;  these  modifications  are  described  in  the  FLAP3  User's 
Guide  which  begins  on  page  19 . 

3.6  Potential  Flow  Calculation 

The  fourth  phase  of  the  complete  three-dimensional  code 
package  utilizes  the  FLAP3  subroutine  which  carries  out  a potential 
flow  analysis  for  the  purpose  of  determining  the  mean  flow  field  in 
the  rocket  cavity.  As  in  the  case  of  the  acousto-modal  analysis, 
only  the  smallest  repeating  geometrical  segment  need  be  considered; 
for  this  calculation,  only  the  cavity  portion  of  the  segment  with  the 
proper  boundary  conditions  is  considered.  The  same  general  formulation 
of  the  finite  element  model  equations  of  motion  is  utilized  except 
that,  in  this  case,  the  fluid  is  considered  to  be  Incompressible.  The 
mass  flow  into  the  cavity  from  the  burning  propellant  surface  is 
modeled  as  a cavity-solid  propellant  Interface  nodal  point  quantity. 
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n is  i-.i  i rul.it  t‘d  by  summing  the  interf  surface  .irc.is  .ir..-iOclaLid 
with  c.u  ii  iiod.iL  point  lyiiif;  on  tlie  c.avi  ty-sol  i d propel  l.inl  interfai:e. 
The  solution  of  the  resulting  set  of  line.ar  equ.iLions  for  the  mean 
flow  veh'city  comiionents  (constant  for  each  cavity  ti'traliedron)  is 
explained  in  Reference  4;  it  is  not  considered  nei:ess.iry  to  go  into 
it  in  det.iil  here. 

3. 7 Evaluation  of  Stability  Integrals 

The  final  phase  of  the  code  consists  of  the  calculation  of 
the  stability  integrals  associated  with  the  various  drivlng/datnping 
coupling  mechanisms  which  occur  in  the  cavity  chamber  in  tlie  presence 
of  combustion  and  flow.  The  stability  integrals  presently  incorporated 
into  FLAP3  are  those  derived  [9)  for  the  three-d linens ion.al  case,  along 
with  the  flow-turning  formulation.  The  use  of  a linear  pressure  (and 
therefore  constant  velocity)  tetrahedral  element  to  represent  the 
cavity  region  enables  an  exact  evaluation  of  the  stability  Integrals, 
given  the  acoustic  nodal  point  pressures  and  element  velocities  and 
the  mean  flow  element  velocities  from  the  finite  element  solutions. 

At  present  FLAP3  does  not  contain  a routine  for  evaluating  nozzle 
damping  or  particle  damping. 

The  calculation  of  stability  integrals  is  for  the  purpose  of 
evaluating  the  driving/damping  coefficient,  a,  a fact  well-known  to 
the  combustion  community.  A positive  a indicates  an  unstable  mode  of 
oscillation  while  a negative  a indicates  a stable  mode.  The  net  value 
of  a computed  by  FLAP3  is  a summation  of  the  computed  values  of 

(pressure  coupling),  a (velocity  coupling),  a (flow-turning),  and 
Ugp  (structural  damping).  The  value  of  is  obtained  from  the  com- 

plex elgensolutlon  described  in  an  earlier  section;  the  other  three 
a-values  are  obtained  from  the  evaluation  of  the  stability  integrals. 

It  is  known  that  the  pressure  coupling  mechanism  always  drives  the 
acoustic  oscillations,  that  the  velocity  coupling  mechanism  may  either 
drive  or  damp  the  oscillations,  and  that  the  flow-turning  mechanism 
always  damps  the  oscillations.  The  response  factors  are  input  into 
the  program  as  multiples  of  the  stability  integrals  for  the  calculation 
of  the  a's  obtained  from  the  different  coupling  mechanisms.  The 
propellant  grain-dependent  response  functions  are  obtained  from  other 
analyses  and  utilized  as  direct  input  into  FLAPS. 
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Section  4 

FLESH3  USER'S  GUIDE 


I,  TITLE  CARD  (80A1) 


Columns 

Variable 

Entry 

1-80 

ITITLE 

Title  information 

Key  word  "END"  stops  job 

II 

, CONTROL  CARD  (4I5.F10.  0, 15) 

Columns 

Variable 

Entry 

1-5^ 

IMAX 

6-10^ 

JMAX 

11-15 

NNC 

Total  number  of  part  boundary  curves 

16-20^ 

IBC 

Number  of  node  code  sequences 

21-30 

SCALE 

Scale  factor  to  multiply  coordinates 
(Default  set  to  1) 

31-35 

NLAY 

Number  of  sections  in  z-direction 

111. 

PART  BOUNDARY 

CURVE  CARDS  (2F10.0,2I5,4F10.0)^ 

Columns 

Variable 

Entry 

1-10 

A(l, 

K) 

X coordinate  (x^  of  ellipse) 

11-20 

A(2, 

K) 

y coordinate  (y^  of  ellipse) 

21-25 

NN 

Number  label  of  curve  (for  tabulated 
curves,  last  card  only) 

26-30 

NP 

Blank  for  tabulated  curve 

Angle  subtended  by  generated  ellipse 

(to  nearest  degree) 

31-40 

AA 

Radius  of  ellipse  along  x-axis 

41-50 

BB 

Radius  of  ellipse  along  y-axis 

51-60^ 

THl  1 

Beginning  and  end  angles  in  degrees, 

61-70^ 

TH2  1 

counterclockwise  positive  from  x-axis 
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IV, 

NODE 

CODE  SEQUENCE  CARDS  (515) 

Columns 

Variable 

Entr 

1-5^ 

IC 

Node  sequence  type 

6-10 

11 

11-15 

J1 

16-20 

11 

(or 

12) 

21-25 

J2 

(or 

Jl) 

V. 

PART 

DEFINITION  CARDS  (1415)^ 

(A  value  of  -1  will  end  this  sequence) 


Columns 

Variable 

Entry 

1-5 

NP 

Part  number  (do  not  define  degenerate 
parts) 

6-10 

11 

Minimum  index  of  part 

11-15 

Jl 

Minimum  index  of  part 

16-20 

12 

Maximum  index  of  part 

21-25 

J2 

Maximum  index  of  part 

26-30 

LI 

Number  label  of  curve  from  (11, Jl) 
to  (12, Jl) 

31-35 

L2 

Number  label  of  curve  from  (12, Jl) 
to  (12, J2) 

36-40 

L3 

Number  label  of  curve  from  (12, J2) 
to  (11, J2) 

41-45 

L4 

Number  label  of  curve  from  (I1,J2) 
to  (11, Jl) 

46-50 

MT 

Part  designator 

51-55 

NN 

Number  of  additional  parts  in  (I) 
direction  whose  constant  (I)  curves 
are  incremented  by  1 

56-60 

13 

Increment  in  I (13)  default  set  to  1 

61-65 

J3 

Increment  in  J (J3)  default  set  to  1 

66-70 

IR 

Set  to  1 to  invoke  Laplace  relaxer 
for  this  sequence  of  parts 

15 


Columns 


1-10 
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VI.  SECTION  LOCATION  CARD  (FIO.O) 

(The  axis  origin  should  be  placed  at  the 
nozzle  end  of  the  cavity) 

Variable  Entry 

ZO  Distance  along  x-axls  (cavity 

longitudinal  axis) 


PART  1 (FLUID) 

II,  J1  TO  15,  J1  (CURVE  1) 

15,  J1  TO  15,  J5  (CURVE  2) 

15,  J5TO  II,  J5  (CURVE  3) 
l1,J5TOI1,  J1  (CURVE  4) 

PART  2 (FLUID  INTERFACE) 
l1,J5TOI5,J5(CURVE  3) 

15,  J5  TO  15,  J6  (CURVE  2) 

15,  J6  TO  n,J6  (CURVE  5) 

II,  J6TO  II,  J5  (CURVE  4) 

PART  3 (SOLID) 

II,  J6  TO  15,  J6  (CURVE  5) 
I5,J6TOI5,J10  (CURVES) 
15,  J10  TO  II , J10  (CURVE  7) 
II,  JIO  Toil,  J6  (CURVES) 


Figure  4.  Finite  element  mesh  generation 
for  cavity-solid  propellant  by  FLESHI. 


FLESH3  Notes 


1.  Each  node  will  be  identified  in  a sequence  as  this  input  labels 
the  node  type.  No  node  may  be  contained  in  more  than  one  node 
sequence. 

2.  Care  is  necessary  in  defining  part  boundary  curves  to  insure  that 
all  parts  are  bound  by  intersecting  curves  and  that  ellipses  are 
cut  by  other  curves  passing  through  them.  All  part  boundary 
curves,  regardless  of  section  location,  are  labeled  and  defined 
in  Cards  III. 

3.  For  tabulated  curves,  THl  = 0.;  TH2  is  blank. 

4.  Enter  for  IC:  1 if  fluid  node  sequence;  2 if  interface  node 

sequence;  3 if  solid  node  sequence. 
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Each  section  of  a common  z-dimension  is  comprised  of  four-sided 
parts  which  are  identified  by  CCW  listing  of  their  part  boundary 
curves;  one  card  for  each  part  (i'igure  4).  MT  is  the  part 
designator  and  identifies  all  elements  within  that  part  as  being 
of  the  same  material.  Enter  for  MT:  1 if  fluid  part,  2 if 

interface  part,  3 if  solid  part.  Interface  parts  are  those  whose 
elements  all  have  at  least  one  side  on  the  fluid-solid  interface. 
Elements  with  only  one  edge  on  the  Interface  are  not  interface 
elements.  Enter  -1  for  NP  to  terminate  each  set  of  part  defini- 
tion cards  which  define  a section  perpendicular  to  the  z-axis 
(cavity  longitudinal  axis). 

Card  VI  is  placed  at  the  end  of  a series  of  Part  Definition  cards 
which  describe  a section.  There  should  be  as  many  Cards  VI  as 
there  are  sections. 


Section  5 

FLSH3P  USER'S  GUIDE 

The  use  of  FLESH3  is  limited  to  continua  with  consistent  part 
description  from  section  to  section.  Consequently,  when  abrupt  changes 
in  cross-section  configuration  occur,  it  is  necessary  to  run  the 
continua  through  FLESH3  in  segments  and  to  post-process  the  results  to 
combine  the  segments.  A post-processor,  FLSH3P,  is  provided  for  such 
cases.  Minimal  FLESH3  output  manipulation  is  necessary  to  enable  FLSH3P 
to  renumber  all  nodes,  elements,  and  interface  surfaces,  and  to 
redesignate  nodes  which  may  fall  on  the  interface  surface  after  con- 
solidation of  the  segments.  FLSH3P  will  not  generate  interface  sur- 
faces perpendicular  to  the  z-axis  (cavity  longitudinal  axis). 

I.  CONSOLIDATION  CONTROL  CARD  (315)^ 


Columns 

Variable 

Entry 

1-5^ 

ND 

Number  of 

segments  to  be  consolidated 

6-10 

IM 

Maximum  I 

index  of  sections 

11-15 

JM 

Maximum  J 

index  of  sections 

II. 

SEGMENT  CONTROL  CARD  (515) 
("ND"  number  required) 

Columns 

Variable 

Entry 

1-5 

NC 

Number  of 

nodes  in  segment 

6-10 

NF 

Number  of 
segment 

fluid  brick  elements 

in 

11-15 

NS 

Number  of 
segment 

solid  brick  elements 

in 

16-20 

NI 

Number  of 
segment 

quadrilateral  Interfaces  in 

21-25 

NZ 

Number  of 
(NLAY  for 

sections  in  previous 
previous  segment) 

segment 

FLSH3P  Notes 


1.  The  consolidation  control  card  is  input  as  the  first  record  with 
segment  control  cards  and  FLESH3  data  following.  The  segment 
control  cards  for  each  segment  precede  the  segment  data  cards. 

2.  A value  of  0 stops  the  job. 
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Section  6 

FLAPS  USER'S  GUIDE 
I.  TITLE  CARD  (80A1) 


Columns 

Variable 

Entry 

1-80 

ITITLE 

Title  information 

Key  word  "END"  stops  job 

II. 

CONTROL  CARD  (315, FI 0.0) 

Columns 

Variable 

Entry 

1-5 

NOD 

Total  number  of  nodes 

6-10^ 

NTS 

Number  of  interface  surfaces 

11-15^ 

NS  LOT 

Number  of  radial  slots 

16-25 

ANGLE 

Angle  (in  degrees)  subtended  by  smallest 
repeating  segment;  must  be  whole  number 
of  degrees  if  NSLOT  = 0 

III.  FLUID 

PROPERTIES  CARD  (I5,2F10.0) 

Columns 

Variable 

Entry 

1-5^ 

NFB 

Number  of  fluid  bricks 

6-15 

BULK 

Fluid  bulk  modulus 

16-25 

FDEN 

Fluid  mass  density 

IV.  SOLID 

PROPERTIES  CARD  (I5.4F10.0) 

Columns 

Variable 

Entry 

1-5^ 

NSB 

Number  of  solid  bricks 

6-15 

SDEN 

Solid  mass  density 

16-25 

PRT 

Poisson's  ratio 

26-35 

GR 

Storage  modulus 

36-45 

GI 

Loss  modulus 
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V.  MODAL  ANALYSIS 

CARD  - ZERO  HARMONIC  (315) 

Coluinns 

Variable 

Entry 

1-5^ 

ZHAR 

Enter  1 if  zero  harmonic  is  desired; 
enter  0 if  not 

6-10 

IZM 

Enter  number  of  lowest  zero  harmonic 
mode  desired 

11-15 

LZM 

Enter  number  of  highest  zero  harmonic 
mode  desired 

VI.  MODAL  ANALYSIS  CARD  - M/2  HARMONIC  (315) 

• 

Columns 

Variable 

Entry 

1-5^ 

MHAR 

Enter  1 if  M/2  harmonic  is  desired; 
enter  0 if  not 

• 

6-10 

IMM 

Enter  number  of  lowest  M/2  harmonic 
mode  desired 

11-15 

LMM 

Enter  number  of  highest  M/2  harmonic 
mode  desired 

VII.  MODAL  ANALYSIS  CARD  -K  HARMONICS  (315) 


Columns 

Variable 

Entry 

1-5^ 

KHAR 

Enter  number  of 
enter  0 if  no  K 

K harmonics  desired; 
harmonics  are  desired 

6-10 

IKM 

Enter  number 
mode  desired 

of 

lowest  K harmonic 

11-15 

LKM 

Enter  number 
mode  desired 

of 

highest  K harmonic 

VIII.  K HARMONICS  CARD  (815) 

(Omit  this  card  if  NSLOT  = 0 or  if  KHAR  = 0) 


Columns 

Variable 

Entry 

1-5 

KHD 

Enter  1 if 

II 

is  desired;  enter 

0 if  not 

6-10 

Enter  2 if 

K = 2 

is  desired;  enter 

0 if  not 

. (continue  through  all  possible 

K harmonics)^ 

36-40 

IX.  COMBUSTION  PARAMETERS  CARD  (2F10.0) 


Columns 

Variable 

Entry 

1-10 

GSPEED 

Speed  of  gas  leaving  the  burning 
surface 

11-20 

CCM 

Mass  fraction  of  particulate  material 

X.  RESPONSE 

FACTOR  CARD  (4F10.0) 

Columns 

Variable 

Entry 

1-10 

RPC 

Pressure  coupling  response  factor 

11-20 

RVC 

Velocity  coupling  response  factor 

21-30 

RET 

Flow  turning  response  factor 

31-40 

RSD 

Structural  damping  response  factor 

Columns 


XI.  NODAL  POINT  DATA  CARDS  (15 , 3F10. 0, 15) 
axis  origin  should  be  placed  at  the  nozzle  end  of  the  cavity  ) 
Variable  Entry 


Nodal  point  number 


X coordinate 


16-25 


y coordinate 


26-35 


z coordinate 


36-40 


NTYPE 


Nodal  point  identifier 


XII.  FLUID  ELEMENT  CARDS  (1015) 


Columns 


Variable 


Brick  element  number 


6-10 

11-15 

16-20 

21-25 

26-30 

31-35 

36-40 

41-45 

46-50 


Global  node 
point  numbers 
corresponding 
to  element  nodes 


Enter  2 for  an  element  on  the  fluid- 
solid  interface;  enter  1 otherwise 


SOLID  ELEMENT  CARDS  (1015) 
(Omit  these  cards  if  NSB  = 0) 


Columns 


Variable 


Brick  element  number 


6-10 

11-15 

16-20 

21-25 

26-30 

31-35 

36-40 

41-45 


Global  node 
point  numbers 
corresponding 
to  element  nodes 


46-50 


Enter  3 
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XIII.  QUADRILATERAL  INTERFACE  SURFACE  CARDS  (615) 

Columns  Variable  Entry 

Quadrilateral  surface  number 

1 
2 
3 
A 

Enter  1 for  a surface  whose  outward 
normal  Is  positive;  enter  2 for  a 
surface  whose  outward  normal  is 
negative 


1-5 

SNUM 

6-10  ] 

Global  node 

11-15  1 

point  numbers 

16-20  1 

corresponding 

21-25  1 

1 to  element  nodes 

26-30^° 

FLAP3  Notes 

1.  An  interface  surface  is  the  quadrilateral  face  of  a brick 
(quasi-hexahedral)  which  lies  on  the  cavity-solid  propellant 
boundary. 

2.  The  program  is  not  limited  to  the  analysis  of  cavities  having 
narrow  radial  slots.  It  is  applicable  to  any  similar  type 
geometry  such  as  a star  or  cloverleaf  pattern,  or  to  an  annular 
cavity  where  NSLOT  = 0. 

3.  This  is  the  number  of  bricks  (quasi-hexahedral)  making  up  the 
cavity  portion  of  a repeating  segment. 

A.  This  is  the  number  of  bricks  (quasi-hexahedral)  making  up  the 
solid  propellant  portion  of  a repeating  segment.  When  the 
"cavity  only"  option  is  exercised,  NSB  is  input  as  0. 

5.  There  are  three  distinct  types  of  harmonics:  zero,  K,  and  M/2. 

The  zero  harmonic  always  exists;  the  K harmonics  are  given  by 


K = 1,  ...,  J 


where 

J = (for  odd  M) 

M - 2 ,, 

J = — - — (for  even  M) 
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where  M is  the  number  of  radial  slots  or  lobes;  the  M/2  harmonic 
exists  only  when  M is  an  even  number.  When  NSLOT  = 0,  only  the 
zero  harmonic  exists. 


6.  Ail  of  the  remaining  data  (Cards  XI  through  XIII)  can  be  automati- 
cally generated  by  the  program  FLESH3  which  is  specially  developed 
for  use  with  FLAP3. 


7.  Each  nodal  point  lying  inside  the  cavity  region  is  identified  by 
1;  each  nodal  point  lying  on  the  cavity-solid  propellant  Inter- 
face is  identified  by  2;  each  nodal  point  lying  in  the  solid 
propellant  region  is  identified  by  3. 

8.  Element  node  numbers: 


When  the  outward  normal  to  a quadrilateral  interface  surface  is 
positive,  the  surface  is  designated  1;  when  the  outward  normal 
is  negative,  the  surface  is  designated  2. 


\ 


\ \ \ 


CAVITY  , 
REGION  ^ 

W 


/ PROPELLANT  / 
/ REGION 

/ 

/ 

POSITIVE 
OUTWARD 
NORMAL 


QUADRILATERAL  INTERFACE 
SURFACE 


10.  If  FI.ESH3  is  used  to  generate  the  finite  element  mesh  data  for 

FLAP3,  all  quadrilateral  interface  surfaces  will  be  automatically 
designated  as  1 (positive  outward  normal).  In  cases  where  the 
outward  normal  is  actually  negative,  the  designation  change  to  2 
must  be  made  manually. 
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Section  7 
CONCLUSIONS 

Simply  stated,  FLAP3  performs  a ] Inear  acousto-modal  analyslr 
of  the  irrotatlonal  motions  of  an  inviscid,  compressible  fluid  coupled 
to  the  irrotatlonal  motions  of  a nearly  incompressible,  I inearly  isotropic 
viscoelastic  solid  and  a linear  potential  flow  analysis  of  the 
irrotatlonal  motions  of  an  inviscid,  incompressible  fluid,  and  then 
determines  the  effect  of  the  flow  field  and  combustion  on  the  calculated 
acoustic  oscillations.  There  are  obvious  limitations  attached  to  any 
code  which  is  as  basic  as  the  previously  listed  restrictions  dictate, 
but  it  is  felt  that  the  developed  code  presented  here  is  probably  as 
sophisticated  as  the  present  state-of-the-art  warrants.  It  is  viewed 
by  the  developer  as  having  much  potential  as  both  a design  and  a 
research  tool.  As  the  state-of-the-art  in  combustion  technology  advances, 
it  is  felt  that  the  code  can  be  revised  relatively  easily  and  updated 
to  include  the  new  technology;  at  least  it  was  designed  with  that  in 
mind. 


One  of  the  most  attractive  features  of  FLAP3  is  the  ease  of  use, 
which,  hopefully,  is  demonstrated  in  the  examples  found  in  the  appendices. 
Other  extremely  important  attributes  are  the  fact  that  it  is  three- 
dimensional,  that  it  performs  a coupled  cavity-solid  propellant  analysis 
(or,  alternatively,  a cavity  only  analysis)  and  that  all  analyses  are 
contained  in  a single  program.  Features  of  the  code  which  do  not  enhance 
its  reputation  also  exist;  they  too  should  be  pointed  out.  It  is  a 
large  program  requiring  a large  amount  of  storage  and  it  may  require 
long  run  times,  as  is  the  case  with  any  three-dimensional  finite  element 
program.  Presently,  the  entire  program  is  in-core  computation,  but  this 
will  probably  be  modified.  In  certain  Instances,  a two-dimensional 
(axisymmetric)  uncoupled  analysis  provides  sufficient  accuracy;  for 
such  cases,  use  of  the  three-dimensional  code  might  not  have  merit. 

It  is  felt  that  the  demonstrated  attributes  of  FLAP3  far  outweigh 
any  foreseen  disadvantages  and  that  it  can  provide  the  means  of  per- 
forming important  analyses  which  were  impossible  before  its 
implementation. 
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Appendix  A.  DIMENSIONS  OF  VARIABLES 


MAIM  PROGRAM  - FLAP 3 

1)  Ttie  following  variables  have  the  fixed  dimensions  indicated: 
K(16,16),  SUM(l6,l(i),  TR(16,16),  FF(10,10),  SB(10,16), 

ART(12),  N(16),  NHARM(16),  ITITLL(80). 

2)  Tlie  following  variables  have  dimensions  at  least  equal  to  N'OU, 

where  MOD  is  the  total  numlier  of  nodes  whose  coordinates  are 
input:  X,  Y,  Z,  MP , MR,  MT , MZ , NH , SP , MTYPE , SS,  AMR, 

AMT,  ANZ. 

3)  Tlie  following  variables  have  dimensions  at  least  equal  to 

MIS  where  MTS  is  the  number  of  quadrilateral  interface 
surfaces:  IMS,  JMS,  LMS , MMS , ISCODE,  SAREA,  XCOM , YCOM,  ZOOM. 

A)  Tlie  following  variables  have  dimensions  at  least  equal  to 
2*NIS:  1ST,  JST,  KST,  IDT,  IDS,  SVX,  SVY,  SVZ,  STAR. 

5)  The  variables  INT,  JNT , LMT,  and  MNT  have  dimensions  at  least 
equal  to  5*NE,  where  ME  = MFB  + MSB,  where  NFB  is  the  number 
of  fluid  bricks  and  MSB  is  the  number  of  solid  bricks. 

6)  The  variables  PR  and  PL  have  dimensions  at  least  equal  to 
(S,P)  where  S is  one-half  of  the  number  cf  repeating  segments 
and  P is  the  number  of  fluid  nodal  points  plus  tlie  number  of 
interface  nodal  points. 

7)  The  variables  MASS,  STIFF,  and  STM  have  dimensions  at  least 
equal  to  (NUDF,MUDF)  where  NUDF  = (Number  of  Fluid  Nodal 
Points)  + 5*(Number  of  Interface  Nodal  Points)  + 4*(Mumber 
of  Solid  Nodal  Points),  or  NUDF  = P,  when  NSB  = 0. 

8)  The  variable  T has  dimensions  of  (2*NPD,NPD)  where 

NPD  = NUDF  - P.  If  NSB  = 0,  the  dimensions  of  T can  be 
(2*P,P).  If  NSB  = 0 and  KHAR  = 0,  tlie  dimensions  of  T can 
be  (1,1). 

9)  The  variable  KHD  has  the  dimension  of  KHAR.  The  variables 
AL,  BW  and  EIG  have  the  dimension  of  LZM,  LMM  or  LKM, 
whichever  is  larger. 

SUBROUTINES  - FXAP3 

Table  A denotes  the  correspondence  between  subroutine  variables  and 
variables  in  the  main  program.  All  subroutine  variables  must  have  the 
same  dimensions  as  the  corresponding  variables  of  the  main  program. 
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T,\BU;  A.  C0RW;S1'0NUENCE  OF  VARIAlUJiS  - FLAl'3 


Subroutine  Variables 

Main  Program  Variables 

Subroutine  MULTCF 

BUM 

MASS  or  STIFF 

Subroutine  MULTOS 

BUM 

MASS  or  STIFF 

Subroutine  MULT 

DUM 

MASS  or  STIFF 

Subroutine  DF.TF.RM 

S 

S'TM 

Subroutine  EQSOLV 

S 

STM 

FI.SH3P 


1)  'Hie  variables  NP  and  NM  have  dimensions  at  least  equal  to  ISNC 
where  ISNC  is  the  sum  of  the  nodal  points  in  all  of  the 
segments  which  are  being  combined. 


2)  The  variable  C has  dimensions  at  least  equal  to  (3, ISNC). 


3) 


The  variable  MF  has  dimensions  at  least  equal  to  (lO.lSNF) 
where  ISNF  is  the  total  number  of  fluid  bricks. 


A) 


The  variable  MS  has  dimensions  at  least  equal  to  (lO.ISNS) 
where  ISNS  is  the  total  number  of  solid  bricks. 


5)  The  variable  MI  has  dimensions  at  least  equal  to  (6,ISNI) 
where  ISNI  is  the  total  number  of  quadrilateral  interface 
surfaces . 
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Appendix  B.  EXAMPLE  ANALYSIS  OF  A TYPICAL 
BEAM  RIDER  MOTOR  CONFIGURATION 


The  beam  rider  motor  configuration  shown  in  Figure  B has  an  annular 
cavity  containing  a solid  propellant  core,  with  solid  propellant  also  at 
the  head  end.  Because  the  cavity-propellant  geometry  undergoes  abrupt 
changes  from  end  to  end,  the  finite  element  input  data  to  FLAP 3 were 
generated  using  both  FLESH3  and  FLSH3P.  Because  FLESH3  does  not  auto- 
matically generate  cavity-propellant  interface  surfaces  which  are  per- 
pendicular to  the  longitudinal  axis,  these  surfaces  and  those  interface 
surfaces  having  negative  outward  normals  were  manually  added  to  those 
interface  surfaces  which  were  automatically  generated.  It  was  also 
necessary  to  change  the  designation  of  some  of  the  nodal  points  and  bricks 
generated  by  FLSH3P,  for  this  configuration  for  the  preceding  reasons.  The 
changes  are  systematic  and  easily  made  and  are  noted  in  the  Input  to  FLAP3. 

Because  the  cavity  is  annular  (no  radial  slots) , the  angle  subtended 
by  the  smallest  repeating  segment  was  arbitrarily  taken  as  9°.  The  fact 
that  only  the  zero  harmonic  exists  for  the  case  of  an  annular  cavity  was 
recognized,  as  reflected  in  the  input  data.  The  pages  following  Figure 
B contain  reduced  computer  output  for  the  beam  rider  analysis. 
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Appendix  C.  EXAMPLE  ANALYSIS  OF  A TYPICAL  FOUR -SLOT 
FINOCYL  MOTOR  CONFIGURATION 


A typical  four-slot  finocyl  rocket  is  sliown  in  Figure  C.  liecause 
the  cavity-propellant  geometry  undergoes  rather  abrupt  changes  1 rom 
end  to  end,  the  finite  element  input  data  to  KLAJ’J  were  generated 
using  both  KLKSH3  and  FLSH3P.  The  pages  following  Figure  C contain 
tl\e  reduced  computer  output  sheets  for  the  finocyl  analysis. 
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CAKO  I.i  I ITLE. 


FLESH3  - FlNOCrt  CI> 

CARO  11,  COMF  R3L. 

I>nx  UMAX  NNC  I0C  scale  NLAT 

S 3 IC  F 1.  3 

CARO  III.  PART  dOUXOART  CURVE. 


X - COORD  r 
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NN 

NP 

A X-AXIS 

B Y-AXIS 

TMETAI 

THETAP 

3.Q  13 

C.  0 00 

0 

0 

0. 00333 

0.03000 

0.00000 
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3.030 

0.  QDO 

1 

0 
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0.  00  000 

0. 00000 
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0.000 
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0 

0 

0.  croio 
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0. 00000 

0. 00009 

i.on 

3.C0C 

0 
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J.Qlu 
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3 

90 
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0.000 
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90 
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90 
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2.  14T30 
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6 

90 
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2.  34700 
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90.  03033 

0.030 
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90 
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0. 00000 
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0 

0 
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11 

0 
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0 

0 
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12 

0 
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J1 
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1 

1 
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3 

1 
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2 

3 

1 

2 

1 

3 

1 

2 

3 

2 

3 

3 

2 

4 

1 

4 

2 

3 

4 

3 

5 

3 

3 

5 
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2 
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- 1 -Q 


C«RO  41.  SECT 13N  wJ^ATlON. 


Z • 
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0.030 

M 

< 

y 

< 

1 

.lOOOu 
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0. 0000  0 

2 

.93900 

0.30000 

0 . 0 CO  0 0 

3 

1 .18900 

0.00000 

0.00000 

W 

? .14Z00 

0.30000 

0 . 0000  3 

5 

2.34700 

6.  008 lib 

0.  OOtiO  0 

b 

.39864 

.31640 

0.00300 

7 

.92625 
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1. 0003  3 

8 

1.17280 

. 19o00 

0.00030 

9 

».  13807 

.19500 

0.00000 

10 

2.33882 

.195U0 

C.COjQO 

11 

. 0 70  71 

.37071 

0 . 0030  0 

12 

.00397 
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0. coco  0 

13 
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14 
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18 
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CARJ  V 


PA?T  OtFINITION 


■JP  U Jl  12  J2  LI  L>  L<»  IT  Nn  13  J3 


:*Ro  VI.  sEcriJN  ldcaiion. 
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V 
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MI 
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1 

17 

.9390C 

0.  lOOOC 
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1 

18 
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8. 13u0  0 

1 
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8. 1300  0 

2 

20 
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8.13000 

3 

21 

. J 9964 
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1 

22 

.92625 

.15400 
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1 

23 
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2 
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2.13807 

. 195GU 
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2 

25 

2.42508 
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b . 1 30  3 0 

3 

26 

.07071 
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8.13000 

1 

27 

.66397 
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1 

28 

.34075 
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2 

29 
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3 

30 
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3 

C*«0  n.  SECT  ION  location. 
Z - COORD 
1 6 • 26  C 


T 
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MT 
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32 
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33 

1 . 189  00 
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3l 
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36 
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36 
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37 

. 9 26  2t> 

. 1 ShOO 
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38 
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2 
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2 

40 
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3 

41 
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1 

42 
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1 

43 
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2 

4l 
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3 

46 
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1 
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22 

21 

t 

2 
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17  18 

23 

22 
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8 
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28 

3 

13 
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24  25 

30 

29 

3 

14 
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25  2 4 

34  35 

40 

39 

3 
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43 

8 
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1 
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CARO  III.  PA«r  ^OU^^OAR*  COKtTfc. 
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SN 

NP 
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0.  OCO 

0 

0 

0. 80033 

0.  00000 

0. 00000 

0. 00003 

3.030 
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90 
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7 
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0 
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-0 
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2 
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3 
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-0 
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-a 
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-0 
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-a 
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0 . 000  C 0 

19.  9100  3 

3 

21 

.09864 

.01540 
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PL£SM3  - ciNOCrL  ciin 

CARD  II.  CONFRJL. 

IH4X  JHAX  HNC  IbC  SCALE  NLAr 

5 3 10  5 1.  5 

CARD  HI.  =>ARr  3DUX0AR*  CDRVt. 
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b 
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S 
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9 
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10 
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11 
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12 
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-0 
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18 
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19 
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2u 
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21 
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22 
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23 
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24 
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25 

2.43913 
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27 
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28 
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29 
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30 
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CARJ  t.  PA?T  OEPIMIION. 

II  J1  12  J2  LI  L;  (.1  L4  HT  NN  13  JS  H 
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4 2 2 3 3 11  5 2 4 3 -0  -0  -0  -3 

; 3 1 4 2 1 r 12  5 3 -0  -n  -a  -3 

6 3 2 4 3 12  r 2 5 3 -0  -0  -j  -9 

7 4 1 5 2 1 11  12  7 3 -C  -0  -0  -3 

t)  4 2 b 3 12  li  2 7 3 -C  -0  -J  -0 

-1  -0  -&  -0  -3  *0  -3  -B  -0  -0  -0  -0  -0  -3 
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31 
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33 
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3 

It 
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3 
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3 

36 
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37 
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38 
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hO 
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2 

43 

.84075 

. 8 40  7 5 

42. 33500 

3 

44 

1 .5  ISlD 

1.  >1816 

42.33500 

3 

45 

1 .73329 

1.73029 

42.  3350  3 

3 
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3 2 1 3 2 1 5 11  3 -0  -0 

H 2 2 3 3 11  i 2 4 3 -13  -6  -0 

5 3 1 4 2 1 F 12  5 3 -0  -0  -0  -D 
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-.2758J 
- . 2 7 75  6 
-.27629 
-.27512 
-.27449 
.69532 
.69463 
.69325 
. 6 8 97  6 
.69530 
.69462 
.69335 
.68952 
.69508 
.69619 
.69572 
.94781 
.94741 
.94697 
,94782 
.94/51 
.94694 
.94776 

.94777 

.94690 

.24420 

.24425 

.24427 

.24429 

.24443 

.24362 

-.98723 


348.252  HZ 


-.98727 
-.98723 
-.  98732 
-.98719 
-.98736 
-.88221 
-.0824<» 
-.08222 
-.88241 
-.08236 
-.88155 
1.0CL94 
1.00C89 
1.0CC94 
1.00095 
1.80C90 

1 . 0 to  C u to 


400E  NUiBE? 


AoOJSTlC  PRESiU^E 


-.866H4 
-.86660 
-.86  657 
-.86  989 
-.8  6 845 
-.86867 
-.86850 
-.87009 
-.86845 
-.86885 
-.86875 
-.27603 
-.27545 
-.27544 
-.2  7 774 
-.27506 
-.27552 
-.27580 
-.27756 
-.27529 
-.27512 
-.27449 
.69532 
.69463 
.6  9325 
.66976 
.69530 
.69462 
.69335 
.68952 
.695C8 
.69519 
.69572 
.94761 
.94741 
.94697 
.94782 
.94751 


r 


I 


j 

I 


•>3 

56 

5’ 

5d 

61 

62 

66 

67 

71 

72 
7o 
77 
61 
82 
8d 
87 

91 

92 

96 

97 
ICl 
162 
106 
107 
111 
112 
116 
117 


.9<.694 
.9<*77b 
.91.777 
.91.196 
.29426 
.24425 
.24427 
.2h429 
. 2..443 
.24352 
-.96723 
-.98727 
-.98723 
-.98732 
-.96719 
-. 96736 
-.08221 
-.06244 
-.08222 
-.u8241 
-.08236 
-.88155 
1.03C94 
1.00089 
1.00C94 
1.01695 
1.04096 
1.00008 


ZERO  HIRMO'IIC  ^OCE  3 

&LPHPC  = .39 

«4PHOC  = -.33 

AlPHFT  * -3.10 

A4PHS0  = 0.  U 

A..PHA  = -3.14 


t 

I 


harmonic  0 

MCOl  M 

frequency 

ZBi'I.Sll  RA0/3EC 

JEGMINT  1 R 

MODE  NUMBEi 

ACOUSTIC  PRES 

1 

.58778 

2 

.58756 

3 

.58756 

-t 

.58743 

6 

.58778 

T 

.58762 

8 

.58747 

9 

.58757 

It 

• .58777 

12 

.58763 

13 

.58776 

16 

-.00696 

17 

-.0  0 695 

16 

-.00  693 

19 

-.00670 

21 

-.00669 

22 

-.00  695 

23 

-.00671 

24 

-.00  70e 

26 

-.00671 

27 

-.00718 

28 

-.00752 

31 

-.58732 

32 

-.58749 

33 

-.58748 

34 

-.58772 

36 

-.58733 

37 

-.58758 

38 

-.58755 

39 

-.§8759 

41 

-.58735 

42 

-.53742 

43 

-.58714 

46 

-.41116 

47 

-.41182 

48 

-.41335 

51 

-.41116 

52 

-.41183 

53 

-.41360 

56 

-.41135 

57 

-.41140 

58 

-.41257 

61 

.99881 

62 

.99844 

66 

.99683 

67 

.99655 

71 

.99664 

72 

.99847 

76 

-.53617 

llA 


a. 


77 

-.53816 

61 

-.53617 

82 

-.53616 

8-) 

- . 5 3 60 1 

87 

- . 53  90  6 

91 

-.  46121 

92 

-<•.8146 

9o 

-.h8122 

97 

-.48147 

141 

-.48136 

102 

-.48033 

106 

1.00104 

1C7 

1.00194 

111 

1.00104 

112 

1.JS102 

lib 

l.COlOO 

117 

l.OOCOO 

5EGM6NT  1 L 

mode  NU*13E^ 

ACOJSTIC  FRESiJiE 

1 

.58778 

2 

.38756 

3 

.58756 

. 56743 

c 

.58778 

7 

.58762 

6 

.58747 

9 

, 58757 

11 

.58777 

12 

.58783 

13 

.58776 

lo 

-.00  690 

17 

-.  JC695 

Id 

-.00  693 

19 

-.  0 0 670 

21 

-.00689 

22 

-.00  695 

23 

-.00671 

2** 

-.00706 

2b 

-.00671  i 

27 

-.00  716  1 

26 

-.  0 0 35  2 

31 

- 58732 

32 

-.58749 

33 

-.68748 

3^. 

36 

37 

38 
3-J 

41 

42 

43 

46 

47 

48 

51 

52 


•.S6772 

-.58733 

-.58758 

-.58755 

-.58759 

-.56735 

-.66742 

-.58714 

-.41116 

-.41182 

-.41335 

-.41116 

-.41183 


115 


iliMlMial 


S3 

-.61360 

56 

-.61135 

57 

-.61160 

56 

-.61257 

61 

.99881 

62 

.99866 

6o 

.99883 

67 

.99655 

71 

.99866 

72 

.99667 

7o 

-.53817 

77 

-.53816 

ai 

-.53817 

82 

-.53816 

66 

-.53801 

87 

-.53966 

91 

-.68121 

92 

-.68166 

9c 

-.<.8122 

97 

-.68167 

ICl 

-.66136 

102 

-.  <.6033 

10c 

1.00106 

10’ 

1.00C96 

111 

1.00106 

112 

1.0C102 

11c 

1.00100 

117 

l.OOCOO 

Hfi-iMO'tiO  MODE  H 


Al®HPC  = 

1.  )5 

Att-rtVC  = 

-.17 

AlPMPT  = 

- 3.  * 3 

AlPHSO  = 

0.  ] : 

Alpha  = 

-2. .5 

116 


i 


I 


n 


harmonic  0 "Out  6 

'RiQUENCt  i87l.967  RAO/jEC  51o.£38  H7 

SEGHclNI  1 R 


node  N0HBE<  acoustic 


3 


3 

11 

12 
13 
It 
17 
13 
13 
21 
22 
23 
Zf 
26 

27 

28 

31 

32 

33 
3h 

36 

37 

38 
33 

41 

42 

43 

46 

47 

4d 

61 

52 

53 
5t 

57 

58 
61 
62 
66 
67 

71 

72 
7c 


-.23375 
-.23350 
-.2  335 1> 
-.23275 
-.23376 
-.23352 
-.23346 
-.23279 
-.23375 
-.23363 
-.23362 
.10233 
.10253 
.1j2E1 
.10338 
.10233 
.10256 
.10246 
.10  364 
.10227 
.10258 
. 1C260 
.14249 
.14291 
.14353 
.14540 
.14251 
.14297 
.14355 
.14539 
.14262 
. 14254 
.14209 
-.0  8866 
-.08609 
— . 0 87  Oh 
— . 0 8666 
-.08814 
-.0  6 688 
-.08853 
-.08856 
-.0  8 751 
-.32376 
-.32358 
-.32377 
-. 32362 
-.32368 
-.32275 
.67705 


117 


77 

81 

8^ 

8a 

87 

91 

92 
9c 
97 

101 

1C2 

106 

107 

111 

112 

11c 

117 


.67725 
.67705 
.07722 
.67713 
.67635 
-.91628 
-.91619 
-.91628 
-.91619 
-.91633 
-.  91668 
1.00067 
1.0C068 
1.00067 
1.00061 
1.00668 
1 • u 0 0 C C 


SESMfcNT  1 L 


NODE  NUIBER  ACOUSTIC  PRESSJv5 


1 

2 

3 

<• 

6 

7 

8 
9 

11 

12 

13 

16 

17 

16 

19 

21 

22 

23 

2h 

26 

27 

28 

31 

32 

33 
3*. 

36 

37 

38 

39 

41 

42 

43 
46 

- 47 

48 

51 

52 


-.23375 
-.23350 
-.23350 
-.23275 
-.23376 
-.23352 
-.23346 
-.23279 
-.23375 
-.23363 
-.23362 
•1C233 
.10253 
.10251 
.10338 
. 16233 

• lu  256 
.10246 
.10364 
.10227 
. 10258 

• Iw  26u 
.14249 
.14291 
.14353 
.14546 
.14251 
.14297 
.14355 
.14539 
, 14262 
.14254 
.14209 

-.08666 
-.08609 
-.08704 
— . 0 8666 
-.06814 


I 

I 


63 

56 

57 
66 
61 
62 
bo 
67 

71 

72 
7c 
77 
81 
62 
85 
87 

91 

92 
9o 
97 

101 

102 

10c 

107 

111 

112 

116 

117 


2£R0  HAR.M0NIC 


ALPrtPC 

s 

2.31 

AlPHOC 

= 

-.2  5 

ALPHPT 

= 

-4.27 

ALPHSO 

= 

£.;c 

Alpha 

s 

-2.21 

-.08683  J 

-.08853  j 

-.U3B56  i 

-.08751 
-.32376 
-.32366 
-.32377 
-.32362 
-.32388 
-. 32275 
.67706 
.briZ'i 
.67705 
.67722 
.67713 
.67635 
91628 
-.91619 
-.91628 
-.91619 
-.91633 
-.91668 
1.0C157 
1.C0C68 
1.  00067 
1.00061 
1.  00C68 
l.OCOOO 

MODE  5 
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